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■ It is denionstrated that the so-called "unavoidable quantum anomalies" can be avoided 
I in the framework of a special non-linear quantization scheme. In this scheme, the quantized 
• hamiltonians are represented by non-linear but homogeneous operators in Hilbert space. 

0^ I The nonlinear terms are of the same order as quantum anomalies, and their role is to cancel 

■ anomalies. The quantization method proposed is applicable to integrable classical dynamical 
PLhI systems and the result of quantization is again an integrable (but, generally, non- linear) 

"quantum" system. A simple example is discussed in detail. Irrespective of the existence 

Ch I of possible physical applications, the method provides a constructive way for extending the 

^ ■ notion of quantum integrability to non-linear spectral problems and gives a practical tool for 
buiding completely integrable non-linear spectral equations in Hilbert space. 
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1 Introduction 



It is known that quantization does not generally preserve zero commutation relations between 
observables. In other words, classically integrable systems may cease to be integrable after their 
quantization. Of course, the quantization procedure is not unique. There are infinitely many 
ways to quantize a system and, for this reason, the integrability can sometimes be recovered 
after an appropriate choice of a quantization scheme. However, the typical situation (especially 
for multi-dimensional systems) is that there are no quantization schemes at all in which such a 
recovering would be possible. In such a case we speak of the "unavoidable quantum anomalies" . 

In order to understand the reason for the existence of unavoidable quantum anomalies, con- 
sider a certain completely integrable A^-dimensional classical system with N mutually commuting 
integrals of motion Hn, n = 1, . . . ,N. Denote by Hn \ n = 1, . . . , N the quantum versions of 
these classical integrals of motion obtained within some fixed quantization scheme. Assume that 
operators Hn^ do not commute with each other. For a randomly chosen quantization scheme 
this is a typical situation. Computing the commutator of these operators, we obtain 



[Hn \ H^^] — h pjiJ^ -\- h -|- . . . . (1-1) 

Obviously, this expansion should start with Otherwise we would obtain non-zero Poisson 
brackets in the classical limit, which is impossible because H^, n = 1, . . . ,N are, by assumption, 
integrals of motion. Now remember that the chosen quantization scheme is not unique. This is 
a trivial consequence of the fact that, instead of operators Hn \ n = 1, . . . , N, we could consider 
arbitrary operators of the form 

Hn = Hj^'^ +hHi^^ + h'^HP + n=l,...,N. (1.2) 

which, evidently, would have the same classical limit. This fact suggests to improve the com- 
mutation relations ( |l.lD by replacing the operators Hn^ by more general operators ( |l.2D with 



appropriately chosen corrections Hn\Hn\---- Requiring that [Hn,Hm] = 0, using ( PI) and 
(|l.2| ) and collecting the terms of the same order in ^, we obtain the system of recurrence relations 
for corrections of our interest. For example, the equation determining the first correction reads 

+ =0, n>m, n, m = 1, . . . , iV. (1.3) 

However, even this first equation clearly demonstrates that the solution of the problem does not 
generally exist. Indeed, the number of unknowns in equation ( |1.3| ) is N . At the same time, the 
number of equations is A^(A^ — l)/2. If > 3, i.e., for more than three dimensional models, the 
system becomes overdetermined and generally has no solutions. This is what we meant saying 
that for multi-dimensional models the situation with quantum anomalies is typical and cannot 
be avoided in the framework of standard formalizm of quantum mechanics^ 

The aim of the present paper is to demonstrate that construction of quantization schemes 
free of any quantum anomalies (at least on the states) becomes possible if one drops out the 
condition of linearity of quantum mechanics. The resulting non-linear Schrodinger equations are 

^Of course, in the presence of some symmetry in a model, the equations for corrections may cease to be 
independent and their effective number may decrease. In this case the quantization without anomalies may 
become possible. Such situation is usually realized in models obtainable in the framework of the r-matrix method 
if the underlying symmetry is su[2) or sii(3). As much as we know, for higher symmetries the problem with 
quantum anomalies in r-matrix approach remains still unsolved. 
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very simalar to Doebner - Goldin equations recently derived in the framework of a rather 
general quantization scheme based on the use of infinite-dimensional algebras of vector fields 
and group of diffeomorphisms |Q, |3|. Speaking of non-linear Schrodinger equations and non- 
linear quantum mechanics, we mean a hypothetic quantum theory in which the observables 
are represented by non-linear operators in Hilbert space. At the present time there are many 
examples of such theories (see e.g. refs. Q ^, ||] and references therein). Because up to 
now there are no experemental indications that the quantum world is non-linear, the measure of 
non-linear effects in all reasonable theories should be very small. In our scheme this condition is 
automatically satisfied because in this scheme the role of non-linear terms is to cancel quantum 
anomalies, and therefore, the order of these terms is the same as the order of anomalies, i.e., at 
most h'^. 

In this paper we restrict ourselves to discussing a particular (but very important) case of 
quantization of completely integrable classical systems. We leave the consideration of the gen- 
eral case to later publications. In the case of completely integrable systems, the procedure of 
quantization consists of three steps: 1) separation of variables in an integrable classical system^ 
2) standard quantization of one-dimensional equations so obtained, and 3) reconstruction of an 
integrable quantum system from the quantized one-dimensional equations. 

The first step is most non-trivial from the practical point of view, because, in practice, it is 
often very difficult to check whether a given classical system is integrable or not, and even if 
it is, in which canonical coordinates it is separable. However, this difficulty can be avoided by 
using the "classical inverse method of separation of variables" (classical IMSV), when instead 
of checking the integrability or separability of a given system, we are simply building systems 
which are separable and integrable by construction |^ |9| . 

The last step is nothing else than the so-called "quantum inverse method of separation 
of variables" (quantum IMSV), which, up to now, has been applied to systems of separated 
equations with linear dependence on separation constants (see e.g. refs. |Tl|, 12|). For such 



equations the resulting quantum systems are always linear. The non-linear quantum operators 
appear only if the separated equations depend on separation constants non-linearly. 

These two classical and quantum versions of IMSV will be discussed in detail in the two 
subsequent sections § and |^. In section ^ we explain what we mean of non-linear quantization. 
The last section |5| is devoted to discussion of specific two-dimensional models for which all the 
aspects of the proposed quantization scheme become transparent and clear. Unfortunately, these 
models, because of their two-dimensionality, are not the examples of models with true quantum 
anomalies. The discussion of the latter would force us to work in at least four-dimensional space, 
which, of course would lead to very cumbersome expressions, non-desirable in the framework of 
this paper. 



2 Classical version of IMSV 



Let Fi = Fi{x, y,zi, . . . , z^r), i = 1, . . . , N he N (may be partially or completely coinciding) 
polynomial functions of + 2 complex variables x,y,zi,... ,zn. Consider N ordinary first- 

"^The separation of variables is understood here in the generalized (Sklyanin) sense. Following Sklyanin we 
call a A''- dimensional classical system separable if there exist such a canonical transformation to canonically con- 
jugated variables Pi,qi, i = 1, . . . , N in which the Hamilton - Jacobi equation for the system becomes equivalent 
to a system of one-dimensional multi-parameter spectral equations Fi{qi,pi,hi, . . . ,hN), i — 1, . ■ . , N with some 
polynomial functions Fi. For more details see ref. [^. 
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order differential equations 

i^, (g.,^|^,/ii,...,/i7v) =0, i = l,...,N (2.1) 

for N functions Si{qi), i = 1, . . . , N. Obviously, all these equations can be explicitly integrated 
after solving the algebraic equations Fi{x, y,zi, . . . , zn) = with respect to the second variable, 
y = Yi{x, zi^ . . . , zn)- We shall write the result in the following form 

S^{qi)= JY^{q^,hu...,hN)dqi, i=l,...,N (2.2) 

stressing the fact that all solutions depend on N arbitrary numbers hi, . . . ,hN parametrizing 
equations (2.1). 

Let us now note that the equalities (|2.l| ) are not violated if we replace the functions Si{qi) 
in them by their sum 

N 

Siq) = J2Skiqk). (2.3) 

k=l 

This gives us the following new system of equations 

i^. (ft,^|^,/ii,...,/ijv) =0, i = l,...,N, (2.4) 

for a single function S'(q) of N variables qi, . . . ,qj\f. This system admits (by construction) a 
total separation of variables and its solution is given by formulas (^^) and ( |2.2D . 

At the same time, the system ( ^.4| ) can be interpreted as a system of algebraic equations 
for N unknown variables hi, ... , hj^j . Solving it with respect to hi, ... , hj^j we obtain 

h^ = Ho.{ci,^^y a = l,...,N, (2.5) 

where i?Q(q, p) are some functions of two A^-component vector variables q = {qi, . . . ,qN} and 
p = {pi, . . . ,p7v}- Interpreting these variables as canonically conjugated coordinates and mo- 
menta, {qa-iPp} = ^af}-, let us consider a A^-dimensional classical dynamical system with the 
hamiltonian 

H{ci,p)=E{HM,p),...,Hr{ci,p)) (2.6) 

in which E{hi, . . . , /iat) is an arbitrarily fixed function of hi, . . . , h^. 

It is not difficult to see that the stationary Hamilton-Jacobi equation for the hamiltonian 
^(q,p), 

.(.^).., 

is separable by construction and its complete solution (i.e. a solution parametrized by A^ arbi- 
trary parameters) has the form 

N 

E = E{hi,...,hN), S{ci) =Y,I Yi{Quhi,...,hN)dqi. (2.8) 

i=l 
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This immediately follows from the fact that, after fixing the form of the function E{hi, . . . , h^), 
the equation ( p.TD becomes equivalent to the separable system (|2.6D . 

So we have demonstrated that any set of functions Fj(x, y,zi, . . . , zn), i = 1, . . . N gen- 
erates a certain separable classical dynamical system by means of the procedure which we shall 
refer below to as the classical inverse procedure of separation of variables. 

It is not difficult to show that the classical system obtained in such a way is not only sepa- 
rable but also completely integrable in the sense that its hamiltonian admits enough mutually 
commuting integrals of motion. In order to demonstrate this fact, it is sufficient to show that 
functions Ha{c[, p), a = 1, . . . , N form a commutative family. 

Indeed, if Ha = Ha{(i, p),a = 1, . . . N, then we have by definition 

F, = F,{qi,pi,Hi,...,HN)=0, i = l,...,N. (2.9) 

Therefore 

p=p^j:?i^p^o. (2.10) 

dqk OQk fr[ oHi dqk 



dFi dFi A dFi dHi 



A ^ +E777f7r^ = 0' (2.11) 
dPk opk fr{ oHi dpk 



and, consequently. 



N 



dFi dHi OF 



^ dHi dqk dqk 



^ dHi dpk 



m 

dpk 



On the other hand, we have 



= {Fi,,FiJ 



N 

k=l 

N 

E 

k=l 



dqu dpk 
dFi, dF,, 



N 

E 

k,l=l 
N 

E 



dqk dpk 

dF,, dF,, dHi 
dqk dHi dpk 

dF,, dHi dFi 



dpk dqk 

dF,, dFi, 
dpk dqk 



^2 



dHi dqk dpk 

N 



dFi, dFi, dHi 

dpk dHi dqk 

dFi, dHi dFi, 
dHi dpk dqk 



+ 



+ 



+ 



^ dF,, dF, 



kl=l 



dHk dH, 



(2.13) 



(2.14) 
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Applying ( ^l2|) and (|2.13D to (|2l|) we obtain 



iV 



dHk dHi 



dqk dpk dpk dqk 
Since the functions Fi^ and Fi^ depend on different variables if ii 7^ ^2, we have 



N 

E 



OF, dFi 



11 12 



{Hk,Hi} = 0, 



for all ii and i2, and, because of the invertibility of the matrix \\dFi/dHj^\\, we have 

{Hk,Hi} = 0. 



(2.15) 



(2.16) 



(2.17) 



It is not difficult to see that the functions Hi , . . . , i^Af constructed in such a way are functinally 
independent and thus the classical system with hamiltonian H = E{Hi, . . . ,Hn) is completely 
integrable. 



3 Quantum version of IMSV 

Let Fi = Fi{x, y,zi, . . . , z^), i = 1, . . . , N he now (may be partially or completely coinciding) 
polynomial functions of A + 2 generally non-commuting variables x,y, zi, . . . , z^- We assume 
that these variables are distributed in the expression in that order in which they are written. 
Consider A*" linear differential equations 

Fi(^q„ih-^,hi,...,hN^^i{qi) =0, i = l,...,N (3.1) 

for N functions "^i{qi),i = 1, . . . ,N. We shall call equations of the type ( ^.l]) multi-parameter 
spectral equations. The role of the spectral parameters in them is played by the numbers 
hi, . . . ,hi\f. The problem is to find all admissible values of these parameters for which the 



system (3.1) has solutions belonging to a certain a priori given classes Wi of functions ^i{qi) € 
Wi, i = 1, . . . , N. The set of all admissible "A-plets" {hi, . . . , hj\f} will be called the spectrum 
of the system (|3.lD . Below we shall assume that the classes Wi, i = 1, . . . , N are chosen in such 
a way that the system ( p.l[ ) has a discrete spectrum. The corresponding discrete set of solutions 
we represent as 

^i{qi) =Ei{qi,hi, . . . ,hN), i = l,...,N, (3.2) 

stressing their correspondence to the values of spectral parameters hi, . . . , h^- 

Let us now note that the equalities ( p. 2D are not violated if we replace the functions ^i{qi) 
in them by their product 

N 



^(q) = n^«(90 (3-3) 

i=l 

belonging to the class W = ^j^iWi. This gives us the following new system of equations 

I— 
dq 



FWg„i/i^,/ii,...,/ijv)^'(q)=0, i=l,...,N, (3.4) 
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for a single function ^'(q) of the A^-component vector variable q = {qi, . . . , qjv}. This system 
admits (by construction) a total separation of variables and its solution is given by formulas 
and (P). 



At the same time, the system (3.4) can be interpreted as a system of algebraic equa- 
tions for N unknown spectral parameters hi, . . . ,h]\[. Solving it with respect to hi, . . . ,h]sf and 
multiplying the result by ^'(q) we obtain 

ha^{ci) = H^(ci,ih^)^{ci), a = l,...,N, (3.5) 



5q 

where Ha{<i,p) is some formal writing for, generally, non- linear differential operators built from 
two A^-component non-commuting operators q = {qi, . . . , q]\f} and p = {ihd/dqi, . . . , ihd/dq]\f}- 
We interpret these operators as operators of coordinates and momenta satisfying the Heisenberg 
commutation relations [q (8> p] = ihl. 

Consider the A^-dimensional non-linear operator 

H{ci,p)=E{Hi{ci,p),...,Hr{q,p)) (3.6) 

where E{hi, . . . , /ijv) is the same polynomial function of N variables hi, ... , h^ as in the previous 
section. At this point we only note that if the arguments of this function are non-commuting 
variables, they should be distributed in the expression in that order in which they are written. 
The product of non-linear operators will be understood hereafter as their composition. 
It is easy to see that the non-linear spectral equation 

H[ci,ih^^{ci)=E^{c^), (3.7) 

is separable by construction and its complete solution (i.e. a solution parametrized by spectral 
parameters) has the form 

N 

E = E{hi,...,hN), ^(q) = n ^i[a](^"'^i'---'^iv). (3.8) 
This immediately follows from the fact that, after fixing the form of the function E{hi, . . . , hj^), 



the equation ( p.7D becomes equivalent to the separable system (3.4). 

So, we have demonstrated that any A^ arbitrarily chosen functions Fi{x,y,zi, . . . ,zn), r = 
1, . . . N generate a certain, in general, non-linear separable spectral equation by means of the 
procedure which we shall refer to as the quantum inverse procedure of separation of variables. 
The equation (|3.7| ) will be called the (non-linear) Schrodinger equation. 

Now we show that the operators Hq, = Ha{q, p) can (in some sense) be considered as integrals 
of motion of a certain completely integrable non- linear "quantum" system. Indeed, from the 
definition of these operators it immediately follows that they are homogeneous operators of order 
one, i.e., for any constant c 

i7„(c^) = cHa'if. (3.9) 
Moreover, as follows from (|3.5| ), these operators have a common set of eigenf unctions: 

^^^ = /i„^, a = l,...,A^. (3.10) 
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By using these two properties, we can consider two chains of equahties: 



Ho,Hp^ = Ho,hp^ = hpHo,^ = hpK^ (3.11) 

and 

HpH^^ = HpK^ = KHp"^ = hahfs^. (3.12) 

Subtracting ( |3.12| ) from ( 3.11| ), we find that 

[H^,Hf,]^ = 0. (3.13) 

Note that the operators Ha are non-linear only if the multi-parameter spectral equations ( p.2| ) 
depend on their spectral parameters non-linearly. In the case of linear dependence on spectral 
parameters the operators Ha will obviously be linear. Assume that the set of solutions of 



equation (|3.7|) is complete in the sense that W (see above) is a Hilbert space and any function 
from W can be expanded in solutions of equation (|3.7| ). If the operators Ha are linear, then we 
can claim that equality ( p.l3 ) holds for all elements of W and thus the operators Ha commute in 



the strong operator sense. In the case of non-linear operators Ha this reasoning does not work, 
and the commutativity should be understood in the in the weak sense, i.e. on the solutions of 
the spectral problem (^). 

4 The non-linear quantization method 

In the previous section we used the term "Schrodinger equation" and the adjective "quantum" 
in order to stress the fact that the equations (^]^) and ( |3.7|) are related to each other by some 
"quantization procedure" . This means that taking in the "quantum Schrodinger equation" ( f3.7] ) 



the classical limit h ^ 0, we obtain the classical Hamilton-Jacobi equation (2^). In order to 
see this, it is sufficient to represent the "wavefunctions" ^'(q) in the form 

^(q)=exp(^^5(q)), (4.1) 

after which it becomes clear that the factorizability of the wavefunctions implies the decompos- 
ability of their logarithms, i.e. classical actions. This enables one to write 

MQi) = (^^SiiQiij, i = l,...,iV (4.2) 

which reduces the problem of verifying our assertion to checking that in the classical limit the 
multi-parameter spectral equations ( pTl]) reduce to the equations ( |2.1D . But this is obvious 
because the substitution of ( [4. 2]) into (|3.1| ) gives an expression whose constant term in h exactly 
coincides with ( |2.1| ). All other terms of higher orders in h vanish in the limit ?i ^ 0. 

5 An example 

In this section we consider a simple example demonstrating how does the proposed scheme work. 
We start with two F-functions of the form 

Fi{x,y,h,g) =x2 + y2__i_/i2_ ("^3.2 + 1^/, + ^ (5 1) 
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and 

F2{x, y, h, g)=x^ + y^- - (^^x^ + l^h-g. (5.2) 

in which the variables x and y are associated with classical or quantum coordinates and momenta 



and h,g are spectral parameters. We see that, generally, the functions ( |5.1| ) and (5^) depend 
on the spectral parameters non-linearly, but if e = or a = oo, this dependence becomes linear. 
For this reason e and a will play the role of non-linear deformation parameters of a system. Let 



us now construct classical and quantum models associated with functions (5.1) and (|5? 



5.1 Classical case 



According to the results of section ^, functions (5.1) and (5.2) can be used for building the inte- 



grals of motion for a certain two-dimensional completely integrable classical dynamical system. 
Denoting these integrals hy H = H{q, p) and G = G(q, p), where q = {qi, ^2} and p = {pi,P2}, 
and following the general prescriptions of section ^, we can write for them the following two 
elementary solvable equations 

+ " 2^^' " (^'^1 + 1^H + G = (5.3) 

and 

+ - - [^Q2 + l^H-G = 0. (5.4) 

Let us take the first integral H as the hamiltonian of a system. Define for this hamiltonian the 
potential T^(q) = -ff(q, 0). Despite the fact that the hamiltonian cannot be generally represented 
as the sum of kinetic and potential energies, this function still reflects some general spectral 
properties of the model. Indeed, as it follows from equations (^) and (5.3), the hamiltonian 



H = i7(q, p) is an increasing function of p^. Therefore, we have the inequality i7(q, p) > ^(q) 
which means that the volume of the phase space bounded by the level surface -ff (q, p) = h 
should be finite if /i < maxy(q) and infinite if h > maxy(q). 



In the particular case, when the parameter is very large, the system (5.3) - ( ^.4| ) becomes 
linear, and, solving it, we obtain the model with hamiltonian 

H = (5,5) 

which is nothing else than the two-dimensional symmetric harmonic oscillator, i.e., the sum of 
hamiltonians of one-dimensional harmonic oscillators. The second integral of motion, G, reduces 
in this case to the difference of hamiltonians of one-dimensional harmonic oscillators, and thus, 
the commutativity of H and G is obvious. The potential V^(q) = q^/2 of this model does not 
have an upper bound and the motion in this system is always finite. 

If a is finite but e is small, then we still have a linear system of equations for H and G, whose 
solution leads to a little bit more comlicated model with hamiltonian 

H = i4^- (5-6) 

2 a2 + q2 ^ ^ 

Although the commutativity of H with the corresponding G is not so obvious as in the previous 
case, it also can be checked by direct computation of the Poisson bracket. Now the potential 
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of the model is bounded by maxy(q) = and therefore the model (|5.6| ) describes finite 

motion if h < and infinite motion if /i > a^/2. 

In the most general case, when both e and a are finite numbers, the equations ( |5.3| ) - ( ^.4| ) 
become quadratic and have now two different solutions for H. For definiteness we choose that 
solution which is continuously connected with particular solutions ( |5.5D and (^.61). It has the 
form 



Reconstructing the second solution G of this system, and checking its commutativity with H we 
can again make sure that the obtained model is completely integrable. As before, its potential 
is monotonically increasing function of and tends to maxy(q) = a^/2 if q^ — > oo. 

Despite the fact that the model ( ^.T] ) looks rather complicated, its Hamilton - Jacobi equation 
is separable by substitution ^(q) = Si{qi) + S2{q), and reduces to the system of two one- 
dimensional equations for Si{qi) and 5*2 (q'2)- This system is nothing else than the initial system 



p^) - (5.4) with pi = dSi{qi)/dqi, p2 = dS2{q2)/dq2, H = h and G = g. Solving it, we obtain 



the complete integral of the model 



5(q)= j dq,^(^^-?jql + ^h^ + h-g + J dq2^ (^^ - 1^ q^ + + h + g, (5.8) 

parametrized by two arbitrary parameters h and g. From this solution it is clearly seen that 
for h < a^/2 the motion in the system is finite, in full accordance with general reasonings given 
above. 

Concluding the exposition of the classical case, let us stress again the fact that in this case 
there is no principal difference between the models associated with functions ( |5.lD and ( ^.2| ) 
with linear and non-linear dependence on spectral parameters. We have seen that in non-linear 
case the models turned out to be more complicated than in the linear one, but they are still 
ordinary classical models admitting quite standard interpretation in terms of Hamilton - Jacobi 
equation. 

5.2 Quantum case 

Let us now use functions ( |5.1| ) and (^]^) for building quantum versions of the models discussed 
in the previous subsection. We start with multi-parameter spectral equations associated with 
these functions and having the form: 

and 





^ \2 / 2 2 



^ ^ + ^2 " 2^^ " l^'^^ + Ij ^ - *2(<?i) = (5.10) 

It is not difficult to see that the system of these equations will have discrete spectrum if we 
require the square integrability of both functions ^i{qi) and ^2{Q2)- Indeed, let us rewrite these 
equations in the form 



dqi 



+ } = ei*i('7i) (5.11) 
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and 



52 ^ 

-n^-^+uj\H^2iq2)=e2'^2{q2) (5.12) 



with 



2h 

^ = \l^--2 (5-13) 



and 



ei = ^h? + h-g, e2 = + h + g. (5.14) 



We see that formulas ( 5.11| ) and ( ^.121 ) look as ordinary Schrodinger equations for simple har- 
monic oscillators, and therefore, if the functions ^i{qi) and ^2{<l2) are normalizable, have stan- 
dard solutions 



^i(gi) =Wn(V^/%i)exp|^|, ^1(92) =^m(V^/%2)exp (5.15) 
and 

ei=ojh{2n + l), e2=ojh{2m + l), (5.16) 

where by Tin and Tim we denoted the ordinary Hermite polynomials. Comparing formula ( |5.16| ) 
with ( 5.13| ) and (|5.14| ) , we obtain the system of two equations 



J-h" + h-g={2n + l)h^|l-'^, + h + g = {2m + l)h^ 1 - (5.17) 

for h and g. The sum of these equations gives a single equation for h: 



^/i2 + /i=(n + m + l)nyi-^, (5.18) 

which, obviously, has a discrete set of solutions. 

Consider particular cases of equation ( 5.16| ). Let a be large. Then the solution of this 
equation can be written down immediately. It is 

h={n + m+l)n. (5.19) 

We see that the spectrum of the parameter h is unbounded. Let now a be finite and e be small. 
Then we obtain a quadratic equation for h. Solving it and choosing that solution which is a 
continuous deformation of ( 5.19 ), we obtain 



In this case the spectrum of the parameter h is bounded by hmax = a?/2. If the numbers n,m 
are large, then the distances between neighbouring spectral points are small, so that the point 
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hmax = is the accumulation point of the spectrum. The general case of arbitrary e and 
a can be considered analogously. The explicit solution of equation ( 5.18| ) is very complicated, 



because ( 5.18| ) is now fourth order algebraic equation. However, the qualitative behaviour of the 
spectrum is the same as in the last case. As before, the spectrum is bounded by hmax = 
and this is its accumulation point. 

Thus, we have demonstrated the normalizability of functions ^i{qi) and ^2iQ2) as well as 
the discreteness and infiniteness of the spectrum of equations ( [5.9[ ) and ( ^.lOD in the interval 
< h < a^/2. Analogously, it can be shown that for h > a^/2 the functions ^'1(91) and ^2('i'2) 
are non-normalizable and the spectrum of multi-parameter spectral equations (|5.9| ) and ( |5.10| ) 
is continuous, in full accordance with the classical case. 

Let us now try to understand whether the spectral values of the parameter h can be con- 
sidered as the eigenvalues of a certain two-dimensional operator in Hilbert space. In order to 
do this, it is sufficient to follow the general prescriptions of section |3|. Let us multiply the 
first equation ( |5.9[ ) by ^'2(92), and the second on e, (^.lOP , by '^i{qi). Introducing the function 
'l'(q) = ^'i((7i)^2(92) and adding the equations (|5.9D and ( B.lOf ), we get a single equation for 



a* 



2(^q2 + l)/i^vl/(q)=0, 



(5.21) 



containing a single spectral parameter h. Consider first the particular cases of this equation. 

Let a be large. Then the terms proportional to /i^ and /iq^ in ( 5.21| ) vanish and we obtain 
the ordinary Schrodinger equation for a two-dimensional harmonic oscillator 



^{-?i2A + q2}vI/(q) =/iM/(q). 



(5.22) 



Let now a be finite and e small. Then the only term proportional to h"^ in ( |5.21| ) vanishes and we 
obtain again a linear spectral equation but with the weight function q^/a^ + 1. In order to reduce 
the hamiltonian of this equation to the hermitean form, we can redefine the wavefunctions as 
^(q) {c^ /a? + l)"'^/^^'(q), after which the equation becomes 



2 7^2^ 



{-n2A + q2} 



Vq^ + c? 



^'(q) = /i^'(q). 



(5.23) 



Let us now consider the general case of arbitrary a and e. Now the spectral parameter h does 



not enter in equation (5.21) linearly and thus, cannot be considered as an eigenvalue of some 
linear operator. The only thing what we can do in this case, is to solve the equation ( ^.21 ) with 
respect to h, and then, multiplying the result by \I'(q) and exchanging the left- and right-hand 
sides, write the nonlinear differential equation 



:2(q +« 



-1 + 



£2 f-^2^^(q)^q2^(q) 

1 H 

{a? + q2)2^(q) 



1 



^'(q) = /i^(q) 



(5.24) 



Note that the ordinary (linear) quantization of the classical system with hamiltonian H leads 
to the linear Schrodinger equation of the form 



\ 



1 + 



(-n^A + q 
{a? + q2)2 J 



^(q) = /i^(q) 



(5.25) 
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with all necessary permutations of non-commuting operators guaranteeing the hermitian sym- 
metry of the corresponding hamiltonian. Expanding equations ( ^.24| ) and ( ^.251 ) in ti^, it is easy 
to see that the first two terms of these expansions coincide. The difference (the non-linearity) 
appears only in the third term and is of order h'^. This non-linearity is of the type (A^'/^')-^^' 
and is very similar to non-linearities appearing in Doebner - Goldin quantization scheme. 



6 Conclusion 

Of course, the example we discussed here, is rather artifical and hardly has some relation to 
reality (even if one belives that quantum world is non-linear), but it clearly demonstrates the 
idea lying in the ground of our approach. At any rate, irrespective of the physical reasonability 
of models obtainable by means of our non-linear quantization procedure, the method which we 
proposed provides a constructive way for extending the notion of quantum integrability to non- 
linear spectral problems and gives a practical tool for buiding completely integrable non-linear 
spectral equations in Hilbert space. 

An interesting application of this method would be construction of non-linear versions of 
quantum Gaudin models associated with higher Lie algebras. As it was noted in ref. |13|, there 



exist multi-parameter spectral equations whose solutions exactly coincide with Bethe ansatz 
solutions for the Gaudin models. Some of these equations depend on spectral parameters linearly. 
The application of quantum IMSV to such equations leads to ordinary (linear) Gaudin models 
(the simplest sl{2) case of such a transformation was considered in detail in ref. |l^). However, 
if the rank of a Lie algebra is suffuciently high, then, along with the "linear" multi-parameter 
spectral equations, there are equations in which the spectral parameters enter non- linearly. The 
application of quantum IMSV to such equations should lead to completely integrable models 
having the same spectra as the ordinary Gaudin models but realized by non-linear operators in 
Hilbert space. An explicit construction of such models is an interesting mathematical problem 
and we hope to consider it in one of the fortcoming publications. 
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